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Abstract. We are interested in the long time behavior of a two-type density-dependent 
biological population conditioned to non-extinction, in both cases of competition or weak 
cooperation between the two species. This population is described by a stochastic Lotka- 
Volterra system, obtained as limit of renormalized interacting birth and death processes. 
The weak cooperation assumption allows the system not to blow up. We study the exis- 
tence and uniqueness of a quasi-stationary distribution, that is convergence to equilibrium 
conditioned to non extinction. To this aim we generalize in two- dimensions spectral tools 
developed for one-dimensional generalized Feller diffusion processes. The existence proof of 
a quasi-stationary distribution is reduced to the one for a d-dimensional Kolmogorov diffu- 
sion process under a symmetry assumption. The symmetry we need is satisfied under a local 
balance condition relying the ecological rates. A novelty is the outlined relation between 
the uniqueness of the quasi-stationary distribution and the ultracontractivity of the killed 
semi-group. By a comparison between the killing rates for the populations of each type 
and the one of the global population, we show that the quasi-stationary distribution can be 
either supported by individuals of one (the strongest one) type or supported by individuals 
of the two types. We thus highlight two different long time behaviors depending on the 
parameters of the model: either the model exhibits an intermediary time scale for which 
only one type (the dominant trait) is surviving, or there is a positive probability to have 
coexistence of the two species. 

Key words : Stochastic Lotka- Volterra systems, multitype population dynamics, quasi- stationary 
distribution, Yaglom limit, coexistence. 



Our aim in this paper is to study the long time behavior of a two dimensional stochastic Lotka- 
Volterra process Z = {Z}, Zf)t>o, which describes the size of a two-type density dependent 
population. It generalizes the one-dimensional logistic Feller diffusion process introduced in 
[7], [H] and whose long time scales have been studied in details in [3]. 

More precisely, let us consider the coefficients 
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The process Z, called Stochastic Lotka-Volterra process (SLVP), takes its values in (M+)^ 
and is solution of the following stochastic differential system: 



where and are independent standard Brownian motions independent of the initial 
data Zq. The extinction of the population is modelled by the absorbing state (0,0), and the 
mono-type populations by the absorbing sets x {0} and {0} x M^. 

This system p.2p can be obtained as an approximation of a renormalized two-types birth 
and death process in case of large population and small ecological timescale. (The birth and 
death rates are at the same scale than the initial population size). This microscopic point 
of view has been developed in [3] concerning the logistic Feller equation and can be easily 
generalized to multi-type models. The coefficients ri and r2 are the asymptotic growth rates 
of 1-type's and 2-type's populations. The positive coefficients 71 and 72 can be interpreted 
as demographic parameters describing the ecological timescale. The coefficients Cij,i, j = 1,2 
represent the pressure felt by an individual holding type i from an individual with type j. In 
our case, the intra-specific interaction rates cn and C22 are assumed to be negative, modelling 
a logistic intra-specific competition, while the inter-specific interaction rates given by C12 and 
C21 can be positive or negative. In case where C12 > 0, individuals of type 2 have a negative 
infiuence on individuals of type 1, while in case where C12 < 0, they cooperate. Our main 
results proved in this paper require two main assumptions. The first one is a symmetry 
assumption between the coefficients 7^ and c^, 



that we will call "balance condition" ( (j2.5p l. It means that the global rates of influence of 
each species on the other one are equal. In particular, the coefficients C12 and C12 have both 
the same sign. The second main assumption ( (|2.8p ) is required in the cooperative case (when 
C12 > and C21 > 0) and is given by 



which compares the intra-specific to the inter-specific interacting rates. This condition will 
be called the "weak cooperative case" . 

Because of the quadratic drift terms and of the degeneracy of the diffusion terms near 0, 
the SLVP can blow up and its existence has to be carefully studied. We prove the existence 
of solutions to (|1.2p in the competition case and in the weak cooperative case. In the first 
case it results from a comparison argument with independent one-dimensional logistic Feller 
processes. In the general case, the existence of the process {Zt)t and a non blow-up condition 
are less easy to prove. If (|2.5p holds, a change of variable leads us to study a Kolmogorov 
process driven by a Brownian motion and the existence of the process is proved using a well 
chosen Lyapounov funtion. That can be done if (12. 8p is satisfied and we don't know if this 
condition is also necessary to avoid the blow-up. Under conditions (12. 5p and ()2.8p we will 
also prove that (0, 0) is an absorbing point and that the process converges almost surely to 
this point. That means that population goes to extinction with probability one. We will 
show this property in two steps. We will firstly show that the process is attracted by one of 





C1272 = C2171 



C11C22 - C12C21 > 0, 
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the boundaries M+ x {0} or {0} x M_|_. Once the process has attained one of them, it behaves 
as the logistic Feller stochastic differential equation and tends almost surely to (0, 0) in finite 
time. Therefore, our main interest in this paper is to study the long time behavior of the 
process {Zt)t conditioned to non-extinction, either in the competition case, or in the weak 
cooperative case, under the balance condition. 

Let us outline that the long time behavior of the SLVP considerably differs from the one of 
the deterministic Lotka-Volterra system which corresponds to the case where 7i = 72 = 0. 
Indeed, a fine study shows that in this case, (cf. Istas [12]), the point (0,0) is an unstable 
equilibrium and there are three possible non trivial strongly stable equilibria: either the 
remaining population is totally composed of individuals of type 1 (the trait 1 is dominant), 
or a similar situation holds for trait 2 or co-existence of the two types occurs. In particular, 
the population cannot goes to extinction. 

In this paper, we want to describe the asymptotic behavior of the SLVP conditioned to 
non-extinction, thus generalizing the one-dimensional case studied in |3j. This question 
is of great importance in Ecology. Although the limited competition resources entail the 
extinction of the population, the extinction time can be large compared to human timescale 
and certain species may survive for long periods teetering on the brink of extinction before 
dying out. A natural biological question is which type will eventually survive conditionally 
to non-extinction. We will show that in the long-time limit, the two types do not always 
disappear at the same time scale in the population and that a transient mono-type state 
can appear. More precisely, we will give some conditions on parameters ensuring mono-type 
transient states (preserving a dominant trait in a longer time scale ) or coexistence of the two 
traits. The main tool of our study will be spectral theory and the conditions will be obtained 
by comparing the smallest eigenvalues of different killed operators. Indeed, conditioning 
to non-extinction, the process can either stay inside the positive quadrant (coexistence of 
traits) or attain one of the boundary (extinction of the other trait). Let us remark that in a 
work in progress , Champagnat and Diaconis are studying a similar problem for two- types 
birth-and-death processes conditioned to non-extinction. 

The approach we develop is based on the mathematical notion of quasi-stationarity (QSD) 
which has been extensively studied. (See [H] for a regularly updated extensive bibliography, 
[TTt [20] for a description of the biological meaning, [HI [lOl [19] for the Markov chain case and 
[3] for the logistic Feller one-dimensional diffusion). In the latter, the proofs are based on 
spectral theory, and the reference measure is the natural symmetric measure for the killed 
process. We will follow these basic ideas. 

In our two-dimensional SLVP case, the existence of a symmetric measure will be equivalent 
to the balance condition and we are led to study the Kolmogorov equation obtained by 
change of variable. Before studying the conditioning to non-extinction, we will in a first 
step study the long time behavior of the population conditioned to the coexistence of the 
two types (the process stays in the interior of the quadrangle (M^)^ as soon as coexistence 
between the two types holds). Our theoretical results, generalizing what has be done in [3j 
to any dimension are stated in the Appendix. The arguments implying the existence of a 
quasi-stationary distribution are mainly similar. The novelty will concern the uniqueness 
of the quasi-stationary distribution, which is shown to be related to the ultracontractivity 
of the killed process semi-group. We prove that the Kolmogorov process associated with 
the SLVP satisfies this setting and conclude to the existence and uniqueness of the QSD 
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for the Kolmogorov process conditioned to coexistence. To deduce a similar result for the 
process conditioned to non-extinction we need to carefully compare the boundaries hitting 
times, and the extinction time. That will give us our main theorem (Theorem 14. ip on the 
Kolmogorov system. Let us summarize here what does it means coming back to the stochastic 
Lotka-Volterra process. 

Theorem 1.1. 1) Under assumptions 12. 5\) and 12. 8\) . the SLVP is well defined on and 
goes to extinction in finite time with probability one. 

2) The long-time behavior of its law conditioned to non- extinction depends on the starting 
point z and is given as follows. 

• For all z'^ >0, if z = {z^,0), then for all AcR^x M+, 

P(,i,o)(^t G A\To > 0) = (m^ ® 6o)iA), 

where is the unique QSD of the logistic Feller process {Y^^)t defined in \2.1\) . 

• For all ^2 > 0, if z = (0, z^), then for all A C M+ x W^, 

P(o,,2)(Zt G A\Tq > 0) = (do m2)(A), 

where m? is the unique QSD of {Y^)t defined in i2.2\) . 

• There is a unique quasi- stationary distribution m on (]R-|_)^\{(0, 0)}, such that for all 
z = {z\z^) with z^ >0, z'^ > 0, for all A C (M+)2\{(0, 0)}, 

F,{Zt£A\To>0) = m{A), 

where Tq is the extinction time. 

• // Ai, (resp. Xi^i, Xi^2) denotes the positive killing rates of the global population, 
(resp. the population of type 1, of type 2), we get 

— Competition case: Ai > Ai,i + Ai,2 o,nd m is given by 

m = 6o rr? + rxi^ ® 5o- 

Furthermore when > M,i (resp. <), (resp. nn?) is equal to 0. 

In other words, the model exhibits an intermediary time scale when 

only one type (the dominant trait) is surviving. 

— Weak cooperation case: we have two different situations. 

* If Xi > Xi^i for i = 1 or i = 2, the conclusion is the same as in the 
competition case. 

* If Xi < Xi^i for i = 1 and i = 2, then 

m = 5q® + ® 5q + tud, 
where mo is proportional to ui. 

We thus have a positive probability to have coexistence of the two 
species. 

Let us remark that our analysis is not reduced to the 2-dimensional case, as it is made clear 
in the Appendix. All the machinery is still available in any dimension. However the explicit 
conditions on the coefficients are then more difficult to write down. That is why we restrict 
ourselves to the 2-dimensional setting. 



• Jo 

• T2 
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2. Existence of the SLVP and boundary hitting times 

Let us denote D = (R;^)^. Let us remark that dD, x {0} and {0} x are absorbing 
sets for the process {Zt)t, as also {(0,0)}. We introduce 

the first hitting time of {(0, 0)}, 
the first hitting time of M+ x {0}, 
the first hitting time of {0} x M_|_, 
• Tqd: the first hitting time of dD (or the exit time of D). 

Of course, some of these stopping times are comparable. For example 

Tqd <Ti <To ; Tqd < T2 < Tq. 

On the other hand, Ti and T2 are not directly comparable. 

Let us prove the existence of the SLVP in some cases. 

Proposition 2.1. //ci2 > and C21 > 0, then there is no blow-up and the process [Zt)t is 
well defined on M+. In addition, for all x € (M+)^, 

p^(ro < +00) = 1 

and there exists A > such that 

sup E^(e^^°) < +00. 

a;e{M+)2 

Proof. In this competition case, the existence of the SLVP is easy to show, by using a com- 
parison argument (cf. Ikeda-Watanabe [11] Chapter 6 Thm 1.1), and the population process 
{Zt)t does not blow up. Indeed, the coordinate {Z^)t, resp. {Z^)t) can be upper-bounded by 
the solution of the logistic Feller equation 



dYl = yJ-iiY^^dB] + {riYl - cii{Y,^f) dt, (2.1) 

respectively 



dy/ = ^^2Y,^dB^ + (r^Y-^ - C22{Y^f) dt. (2.2) 

These one-dimensional processes have been introduced in [TJ [H] and studied in details in ^ . 
It's easy to deduce (by stochastic domination) that the processes Z^ and become extinct 
in finite time. 

The a.s. finiteness of each Tj, hence of Tgjj, thus follows. It has also been shown in [3] that 
the absorption times from infinity have exponential moments. 



□ 



Let us now consider the general case. We reduce the problem by a change of variables. 



Let us define {Xl,Xf) = (2^^1/71,2^^^/72). We obtain via Ito's formula 



dX^ 



* V 2 8 8 2X1) 



* * V 2 8 8 2X1 
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In all the following, we will focus on the symmetric case where X is a Kolmogorov diffusion, 
that is a Brownian motion with a drift in gradient form as 

dXt = dBt - VV{Xt)dt. (2.4) 

Indeed, all the results of spectral theory obtained in the Appendix have been stated for such 
processes. Obvious computation shows that it requires the following balance condition on 
the coefficients: 

ci2 72 = C2i7i- (2-5) 
This relation is a symmetry assumption between the global interaction rate of type 2 on type 
1 and of type 1 on type 2. (Recall that the coefficients 7^ describe the ecological timescales). If 
(|2.5|) holds, the coefficients C12 and C21 have the same sign, allowing inter-species competition 
(ci2 and C21 > 0) or inter-species cooperation (ci2 and C21 < 0). 

Under this condition, the process X is called the stochastic Lotka-Volterra Kolmogorov pro- 
cess (SLVKP). The potential V is then equal to 

V{x\x') = lyj log(x^) + '-^^^^ - + aix^ix^, (2.6) 

where 



2 ^ V °' ' 16 

i=l,2 



C12 72 C21 71 , , 

a = = . (2.7) 

16 16 ^ ^ 

Let us prove the existence of the SLVKP using the L^-norm as a Lyapunov function, under 
a weak cooperative assumption, that is if 

a < and C11C22 — C12C21 > 0. (2.8) 

Theorem 2.2. Assume balance condition 112. 5\) . weak cooperative assumption Ii2.8\) . then 
there is no blow-up and the processes (Xt), and then [Zt), are well defined on M+. 
In addition, for all x & D, 

^x{TaD < +00) = 1, 

for hath X and Z . 

Hence, under the assumptions of Theorem l2.2l Hypothesis (HI) of Definition lA.3] in Appendix 
is satisfied, what we shall use later. 

Proof. Let us compute 

d\\Xf, = d{Xjf + d{X^f 

= 2{XldB\ + XidBf) + Y^^iXlf (^r, - ^^^^^^ - SlHl^^^ dt. 
The quartic function appearing in the drift term is thus —q{{x^)'^, (x^)^), with 

Decomposing 



{u, v) = cii7i(u)^ -I- C2272('f )^ + 32auv. 



. . , 16a \ . . 

q[u, V) = C1171 u H V H (C11C22 - ci2C2i)7i72 

cii7i J C1171 
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and since a < 0, a necessary and sufficient condition for q{u, v) to be positive on the first 
quadrant (n > 0, f > 0), and to go to infinity at infinity, is tfius C11C22 — C12C21 > 0. Hence 
tlie drift term in the previous S.D.E. is negative at infinity. It easily follows that 

sup E{\\Xtf) < +00, 



ensuring that the processes (Xf), and then (Zf) are well defined on M+. 

Let us now study the hitting time of the boundary dD. We will compare {Xj, X^) with the 
solution of 



dU? = dBj + ( '-^^ - If^ilATl^ _ z±^^Z±^Z±±_ \ dt . 

* * V 2 8 8 y 

Assuming ()2.5p and p.Sp . the diffusion process {U} ,Uf) exists and is unique in the strong 
sense, starting from any point. We consider now the solution built with the same Brownian 
motions as for X. 

We shall see that, starting from the same (x^,x^) in the first quadrant, and for all t < Tqj^, 
X} < Ul and X"^ < U^. 

To this end we can make the following elementary reasoning. Fix lo and some t < Tq£){lo). 
Let us define s Wl = X] — [/] for i = 1, 2 and for s < t. Of course Wq = 0. Due to the 
continuity of the paths, s 1— > Wl is of class and solves an ordinary differential equation 
such that ^TO|^^„ = -^<0. 

Moreover we remark that if at some time u < t, = 0, then ^(W])|^^^ = — < 0- It 
follows that Wg < for < s < t, yielding the desired comparison result. 
Denote Sqd the hitting time of dD for the process U. We thus have Sqd > Tq£). It is thus 
enough to show that Sqd is almost surely finite. But, remark that di^ = e~^^^^dx, with 

Q{x\x^) = Y: ('-^^^^ - ^) + 2a{x')\x'f, (2.10) 

i=l,2 ^ ^ 

is an invariant (actually symmetric) bounded measure for U. The process U is thus positive 
recurrent and it follows that, starting from any point in the first quadrant D, Sqd is a.s. 
finite. □ 



These comparison arguments allow us to obtain other interesting properties of the process 
X, that we collect in the next proposition 

Proposition 2.3. Under the assumptions of Theorem \2.<A the following holds 

• there exists A > such that sup^g^) ¥.x{e^'^^^) < +00, 

• for all X € D, ¥x{Tqo = Ti) > for i = 1,2 (recall that Ti defined at the beginning 
of section 2 is the hitting time of each half axis), and Px{Tqd = Tq) = (recall that 
Tq is the hitting time of the origin). 
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Proof. The first point is an immediate consequence of the same moment controls for both 
Y and U, thanks to the comparison property. We already mentioned this property for Y in 
Proposition 12.11 

The same holds for U since U is known to be ultra-contractive (see definition |B3] in Appendix 
B2) under the hypotheses of theorem 12.21 The ultracontractivity property follows from the 
fact that the invariant measure e~^^^^ dx of the process f/, {Q defined in (I2.10p ). satisfies the 
conditions of Corollary 5.7.12 of |21j . 

One sometimes says that X satisfies the escape condition from D. 

In order to show the second point we shall introduce another process. Namely define for 
i = 1, 2, as the solution of the following stochastic differential equation 

dHl = dBl+('-^- ^^1^^ -^\dt. (2.11) 
* * V 2 8 2Hl) ^ ' 

Of course and iP' are independent processes, defined respectively up to the hitting time 
of the origin. We decide to stick in after it hits 0, as for X*. 

On the canonical space Vtt = C([0, t\,D) we denote by and the laws of the processes 
{Xs/^Tgn)s<t) and {HsATgo)s<t) starting from the same initial point x in D. We claim that 
and are equivalent. This is a consequence of an extended version of Girsanov theory 
as shown in [3] Proposition 2.2. One then have that for any bounded Borel function F defined 
on fit, 

where 

A{t) = a {colfico^f - a{x,f{x2f -a f f 2a [u^lfi^^f {{ulf + {u^^f) - {{u^lf + {u^^f) 



-in + r2) {colfico'sf + \{culi + C2272) (^l)'(^.')' {{u^lf + (u^'sf) + {{coif + {u^f) 

andE^ (resp. E"**-) denotes the expectation w.r.t. toV^ (resp. V^). Remark that ^(t A Tg/j) 
is well defined, so that the previous relation remains true without the ^t<Tgr>{Lo) replacing 
A{t) by A{t A Tqd), i.e. 



E^ \F(uj)] = E^ 



This shows the claimed equivalence. 

It thus remains to prove the second part of the proposition for the independent pair (H^,H^). 
But as shown in [3j Proposition 2.2 again, for each i = 1,2, 

E^' [F(^) = E*^ [F(u;)Ii<ToHe^^*^ 

for some almost surely finite B{t), E^ being the expectation with respect to the Wiener 
measure starting at x'. It follows that the V^^ law of Tq is equivalent to the Lebesgue 
measure on ]0,-|-cx3[, since the same holds for the law of Tq. The law of (Ti,T2) is 
thus equivalent to the Lebesgue measure on ]0, -|-cxD[(g)]0, -|-oo[ yielding the desired result for 
H hence for X. □ 
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3. Existence and Uniqueness of the Quasi-Stationary Distribution for the 

Absorbing Set dD 

We can define a quasi-stationarity notion associated with each absorbing set O, dD, x 
{0}, {0} X M+. The results developed in Section 3 will refer to the absorbing set dD. Its 
complementary in (M+)^ is D, which is an open connected subset of M^, conversely to the 
complementary of other absorbing sets. 

Let us recall what a quasi-stationary distribution is. If F denotes an absorbing set for the 
(M+) ^-valued process Z and Tp the hitting time of this set, a quasi-stationary distribution 
(in short QSD) for Z and for this absorption event is a probability measure u satisfying 

F^{Zt eA\TF>t) = u{A), (3.1) 

for any Borel set A C ]R_|_^\{i<'} and t >0. A specific quasi-stationary distribution is defined, 
if it exists, as the limiting law, as t ^ oo, of Zt conditioned on Tp > t, when starting from a 
fixed population. That is, if for a all x € R_|_^ \ {F}, the limit 

lim F^{Zt e A\TF>t) 

exists and is independant of x, and defines a probability distribution on \ {F}, then it 
is a QSD called quasi-limiting distribution, or (as we will do here) Yaglom limit. 
It is thus well known, (see [8j), that there exists Xp > such that 

p^(rF >t) = e-^^K 

This killing rate A p gives the velocity at which the process issued from the i/-distribution get 
to be absorbed. 

Our aim is now to study the asymptotic behavior of the law of Xt conditioned on not reaching 
the boundary. All the material we need will be developed in Appendix. The later essentially 
extends to higher dimension the corpus of tools introduced in [3j. The spectral theory is 
developed for any Kolmogorov diffusion in M^, and for any dimension d. Then the existence 
of a quasi-stationary distribution is obtained. The novelty is the uniqueness result, since it 
is deduced from the ultracontr activity of the semigroup. 

We will extensively refer to this Appendix, to study the problem of existence of a quasi- 
stationary distribution for the SLVKP, with the potential V defined in (12. 6p . 
We introduce the reference measure, given by 

li{dx\dx^) = e-2^(^')dx = e-«(^''^') dx' dx" 

x^ 

where Q is the symmetric polynomial of degree 4 given in ()2.10p . 

It is the natural measure to deal with, since it makes the transition semi-group symmetric 
in L^(//). One problem we have to face is that the measure has an infinite mass due to 
the behavior of its density near the axes. Remark that the density is integrable far from the 
axes, which is equivalent to 

/ e-^^^'^dx < +00. (3.2) 
Jd 
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This property is immediate in the competition case {ca, 'ji and a are positive) and has already 
been shown in the proof of Theorem 12.21 when (|2.8p holds. 



Proposition 3.1. Assume 112. 5]) and 112. 8\} . Then, there exists some C > such that for all 
X £ D, 

\VV\'^{x) - AV{x) > -C. (3.3) 

We deduce that the semigroup Pt of the SLVKP killed at time Tqo, has a density with respect 
to fi, belonging to Lp{dfi). 

Proof. Under ()2.5p and (j2.8p , we have seen that the explosion time ^ is infinite and thus ()A.3P 
is satisfied. Then the conclusion of Theorem lA.ll holds: one proves by Girsanov's theorem 
that the semigroup Pt of the SLVKP has a density with respect to fi. 
Moreover, computation gives 

|Vyp(x) - AV{x) = 

1 X^ 71(2;^)^ 1, o,o\^ / I X^ 72(3^^)^ / 1n9 g 

^ + riy + cii^\^ + 2axHx'fj + ( " ^ + ^2y + C22^\^ + 2a{x'fx' 
, 1 , n 371 (x^)^ . , 2^2 , 1 , ^2 372(x^)' „ . 1.2 

+2(^ + y-""^^-'"^^) +2(^ + y-^^^^^-'"(^) • 

Hence, we observe that the terms in |Vyp — AV playing a role near infinity are equal to 
f (^ + "?) (when one of the coordinates is close to zero) and to 

{^{x^f + lax.ix^y + (^xi + 2a.2(x^)^ 
= (xr {-'^{x'f + 2a(.^)f + (x^f {-'-^{x^f + 2«(x^)f . 

The two terms in factor of (x^)^ and (x^)^ in the first quantity will not be together equal to 
as soon as a > or as the determinant of the system 

^Yi + 2aY2 = 
2aYi + ^Y2 = 0. 

is non zero, what is satisfied under the condition (j2.8p . It follows that |Vyp(x) — AV{x) 
tends to +00 as |x| tends to infinity. Since |Vyp — AV is a smooth function in D, (j3.3p 
follows. Then, we deduce from Theorem lA.ll that for each t > 0, the density of Pt belongs 
toL2(d^). □ 

Let us now state our first main result. 

Theorem 3.2. Under Assumptions \2. 5\) and \2. <gj] . that is if 

• C1272 = C21 71, 

• ifa<0, C11C22 - C12C21 > 0, 

there exists a unique quasi- stationary distribution ui for the stochastic Lotka-Volterra Kol- 
mogorov process X, which is the quasi-limiting distribution starting from any initial distri- 
bution. 



COMPETITIVE AND WEAK COOPERATIVE STOCHASTIC LOTKA-VOLTERRA 



11 



In particular, there exists Ai > such that for all x ^ D, for all A C D, 

lim e^^'F^iXt G A\T9d > t) = ^A). (3.4) 

t— »oo 

Proof. The proof is deduced from Appendices A and B. Standard results on Dirichlet forms 
(of. Fukushima [9]) allow us to build a self-adjoint semigroup on L^(/x), which coincides with 
Pt for bounded functions belonging to L^(^). Its generator L is non-positive and self-adjoint 
on l?{dfj,), with D{L) D C^iD). Its restriction to C^{D) is equal to 

Lg = ^Ag-V.Vg, geC^{D). 

We now develop a spectral theory for this semigroup (also called Pt), in Lp{dn). 

We check that the hypotheses (H) introduced in Definition IA.3I required to apply Theorem 

IA. 4I are satisfied under the assumptions (j2.5p and (j2.8p . Indeed (H4) is obviously satisfied 
and (HI) and (H2) are deduced from Proposition 13. 1[ Furthermore, using for example polar 
coordinates, one easily shows that Condition (j2.8p implies that 

G{R) = inf{|Vyp(x) - AV{x); \x\ > R and x £ D} > cR^ 

for positive constant c. Thus (H3) holds. 

Since Hypotheses (H) and (HI) are satisfied, Theorem IA.4I implies that —L has a purely 
discrete spectrum of non-negative eigenvalues and the smallest one Ai is positive. Let us 
prove that the associated eigenfunction 7]i belongs to L^{dfi). We have shown in Section B.2 
that r]ie~^ is bounded. Thus 

/ rii{x)dfj.{x) = / r/i(x)e"^^('')(ix < ||7/ie~^||oo / e^^'^^'^dx < +oo, 
Jd Jd Jd 

since e-^W = i e-^«("''"'). 

Thus the eigenfunction -qi belongs to L^((i/i) and therefore, as proved in Theorem \B.2\ the 
probability measure ui = J^^'^^, is the Yaglom limit distribution. 

In order to show the uniqueness of the quasi-stationary distribution, we apply Proposition 

IB. 121 relating this uniqueness property to the ultracontractivity of the semi-group Pt ■ Let us 
show that the sufficient conditions ensuring ultracontractivity stated in Proposition IB. 141 are 
satisfied by the SLVKP. 

The function V is bounded from below in Df, = {y £ D; d{y, dD) > e}. In addition, Condition 
implies that V{R) = sup{y(rc),x G FD, \x\ < R} < c'R'^ and G{R) > cR^, for positive 
constants c' and c. Hence Condition (jB.lOp in Proposition IB. 14) is satisfied with 7^ = k^'^^"^. 
Thus the killed semi-group Pt of the SLVKP is ultracontractive and then the uniqueness of 
the quasi-stationary distribution holds (cf. Proposition IB . 12P . 

Hence existence and uniqueness of the quasi-stationary distribution holds for X. □ 

Remark 3.3. Since the laws of Xt and Zt are related via an elementary change of variables 
formula, a similar result will be true for the stochastic Lotka-Volterra process Z. 

4. Explicit Quasi-stationary Equilibria - Mono-type transient states. 

In the above Section we were concerned by conditioning to coexistence's event. Let us now 
come back to our initial question, that is the long time behavior of the process conditioned 
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to non-extinction. The SLVKP dynamics is particular in the sense that once hitting the 
boundary dD^ the process will no more leave it. Hence, for t > Tq£) the process will stay 
on one half axis, and the dynamics on this axis is given by the process {Hl)t defined in 
()2.1ip . that has been extensively studied in [3]. Thus we know from [3] that for i = 1,2, 
there is a positive killing rate Ai^j > and a unique quasi-stationary measure vi^i on the 
axis x-' = characterized by the ground state r/i^j (eigenfunction related to Ai^j), which is a 
positive function, bounded and square integrable with respect to the corresponding symmetric 
measure /i* on each axis. More precisely, we have 

with 

„,{dx') = e-'ff '^'(^^"V and %(n) = ^ - + 

2u 2 8 

In addition, C Ml, 



lim F^^{Xi £ A\Ti >t) = i^i,iiA), (4.1) 



We deduce from this study that for all x^ > 0, for all A C x (resp. > and 

^{x\o){Xt e A\To >t) = 1/1,1 ® So{A), 
(resp. P(o,x2)(^t G A\To > t) = 6o <S) i^i,2iA)). 

We are now led to study, for x G (M*^)^ and for A C (M^\{(0, 0)}, the asymptotic behavior of 

r (X (^AIT , .X ^AXt^^A) F,{T9n>t) 
P.(X, eA\n>t)- ^^(^r^ P.(ro > t) 

where To is the hitting time of the origin. We have seen in the previous section that under 
(12. 5p and (j2.8p . the hitting time Tqd is almost surely finite and that for any y G dD, Tq is 
¥y almost surely finite too. Let us now study the asymptotic behavior of 

^xjTdD > t) 

Px(To > t) 

We will compare the three different killing rates Ai, Ai,i, Ai,2 corresponding to the stopping 
times Tqu, Ti, T2, (recall that Tj denotes the hitting time of the axis x^ = 0). 

Notice that if ci,2 = 0, (X^X^) = {H^,H'^), Ai = Ai,2 + Ai,i, 1^1 = 1^1,1 (g) 1/1,2- Indeed a 
standard (and elementary) result in QSD theory says that if z/ is a QSD with absorbing set 
C, the Pjy-law of the hitting time Tc of C is an exponential law with parameter A, where A 
is exactly the killing rate. In addition, Fx{Tc > t) behaves like e~^* for large t. Since the 
minimum of two independent exponential random variables with parameters Ai,2 and Ai,i is 
an exponential variable of parameter Ai,2 + Ai,i, we get Ai = Ai,2 + Ai,i. 

The following decomposition is the key point to prove our main theorem. For X G (M+)2\{(0,0)} 
P,(ro >t)= ^^{Tqd >t)+Yl {^dD < t , X^^^ = , P(^^^^^o)(ro > t - Tqd)) (4.3) 

where j = 1 if i = 2 and conversely. We are interested in the asymptotic behavior of this 
quantity, i.e we have to compare the three killing rates Ai, Ai,i and Ai,2- 
We obtain our main result. 
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Theorem 4.1. Under \2. 5]) and \2.^) . there exists a unique probability measure m such that 
for all X £D, for all A C (R+)2\{(0, 0)}, 

lim F,^{Xt £ A\To>t) = iy{A). 

In addition, we have the following description of v. 

• Competition case (ci2 and C21 positive). We have Ai > Ai^i + Ai_2, and the support 
of the QSD in included in the boundaries. 

Furthermore when Ai^2 > Ai^i (resp. <), u is given by vi^i (8 5q (resp. 5q (8) ^^1,2)- 
In other words, the model exhibits an intermediary time scale for which only one 
type (the dominant trait) is surviving. 

• Cooperation case (C12 and C21 negative). We have two different situations. 

— // Ai > Ai_i for i = 1 or i = 2, the conclusion is the same as in the competition 
case. 

— If \i < Xi^i for i = 1 and i = 2, then 
^ J^i,i ®5o + 5o® , , vi,2 + (4.4) 



l + E-,,=iA7|^ VAi,i-Ai-' ^1.2-^- 
where 

We thus have a positive probability to have coexistence of the two species. 

Remark 4.2. The only remaining case is the one where Ai = Ai^i = Ai^2- However the 
proof of the theorem indicates that this situation is similar to the competition case, though 
we have no rigorous proof of it. In the discrete setting (as claimed in [3]), a fine analysis of 
Perron-Frobenius type is in accordance with our guess. 

Proof. Our domination arguments allow us to compare the killing rates : 

• The competition case. This is the case C12 > 0. In this case we can show with a 
similar argument as in the proof of Theorem 12.21 that starting from the same initial 
point, XI < HI for i = 1,2. Hence Ai > Ai^i + Ai^2- Since the killing rates are 
positive, it follows that Ai > Ai^j for i = 1,2. In particular E^; |^e'*'^'*"^'9°] < +00. 
Hence 

lim inf e^i'-* (Tod < t , = , P.^. q) (^o > * - Tqd)) = 

= hminfE, flT,^<tI^, =oe^^-^s^e^i-(*-^«^)P(o^. Jnyt-Tao)) 

> E^ (ix;^ =0 e^i>"'^s^ ViA^ho)) > ^ ( ^* ^^^^t for one i), 

according to Fatou's lemma, the positivity of the ground state and by Proposition 
12. 3i It follows that the rate of decay of Pa;(To > t) is at most e~^i'»*, while the one of 

Px{TdD > t) is e~^^^, hence as t — > +00, 

F,{XteA\To>t)^0, 

if ^ C I?. Hence, the support of the quasi-stationary distribution will be included 
in the boundaries. Thanks to Proposition 12.31 we know that both terms in the sum 
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Y^j^^i 2 in (|4.3p are positive, so that the leading term in the sum will be equivalent to 
^x{To > t). If Ai^i > Ai_2 this leading term is of order e~^^'^'', the proof being exactly 
the same as before. The value of the quasi-stationary distribution follows. 

• The weak cooperative case. This is the case if ci^2 < 0. Here a comparison 
argument gives Ai < Ai^i + Ai^2- But there is no a priori reason for Ai to be smaller 
than Ai^j. In particular if Ai > Ai^j for i = 1 or 2, we are in the same situation as in 
the competition case, and the quasi-limiting distribution is supported by dD because 
one exits from D by hitting = with a positive probability as we mentioned in 
proposition 12.31 

It remains to look at the case Ai < Ai^j for i = 1, 2. 

Denote by ipj the law of Tqd when the process exits D by hitting = 0. Denote 
by Q the conditional law of X^^^ knowing Tqd = s and ^Xgo ~ ^' "^^^^ 

e^i* P, (Teo < t , X^^^ = , P(^^^^ o)(To > t - Tqd)) 

r+oo 

= / e^''ls<tE^^{lTo>t~s)i^j{ds). 
Jo 

It has been proved in Corollary 7.9 of [3J that for any A < Ai^j, 

sup Eg^[e^'^^ < +00 

where 0* describes the set of all probability measures on x-' = 0, > and Tq is the 
first hitting time of for the one dimensional logistic Feller diffusion H^. 

Hence E^i(lTo>t-s) < Ce"^(*-") for some universal constant C and all s. It follows 
that, provided Ai < Ai^j (in which case we choose Ai < A < Ai^j), for all T > 0, 

lim / e^i* I,<t E..(lTo>t-.) i^jids) = . 

It remains to study 

^ hrn^ e^i* {t>T9D>T, = , P(^^^^ ^q)) (^o > t - Tqd)) . 

To this end we first remark that for T large enough we may replace P^; by F,^^ . Indeed, 
denote 

h{T, t,y)=¥y{t-T>T9D>0, X^^^ = , P(^^^^^o)(^o >t-T- Tqd)) ■ 

Since vi is the Yaglom limit related to D and < h{T, t,y) < 1 for all y G and all 
T < t, for e > one can find T large enough such that for all t > T, 

P. {t>T9D>T, = , P(^^^^^o)(^o > t - Tqd)) = 

= E,{lT,^>ThiT,t,XT)) 

K,,{lTso>Th{T,t,XT)) 

= P.1 {t>T9D>T, = , P(^^^^ {To>t- Tqd)) 

where a b means that the ratio a/b satisfies 1 — e < a/b < 1 + e. 
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Now, since i^i is a quasi-stationary distribution, starting from vi, the law of Tqd 
is the exponential law with parameter Ai, so that 'ipj{ds) = Cj Ai e~^^'^ds where cj = 
Pi/j(X;^^^ = 0) is the exit probability first hitting the half-axis j. In addition the 
conditional law Q does not depend on s, we shall denote it by vr* from now on. This 
yields 



e^i* P,, (t>T9D>T, X{^^ = , P(^^^^^o)(To > t - Tqd) 

e^i* / Cj Ai e-^^' P^.(To > t - s) ds 
Jt 

Cj Ai / e-(^i--^^)(*-^) (^e^iAt-s) p^^^j^^ > t - s)'^ ds . 



t 

T 

Recall that z/i j is the unique Yaglom limit on axis i. As shown in [3], for any initial 
law TT, in particular for vr*, lim^^+oo e'*'^''* P7r(To > t) = 1. Using Lebesgue bounded 
convergence theorem we thus obtain that for all T > 0, 



hm c^''F,,(t>T9D>T,X'r, =0,F^^. ,^{To > t - T9D , . , . 

Since the result does not depend upon T, we may use our approximation result of 
Pi by Pi/j for all e hence finally obtain for all x, 

tiToo '''' {t>T3D>T, X^^^ = , P(^. ,0) iTo>t-Ten))= j^^^ • 

Let us now consider A C M^\{(0, 0)}. To compute Fx{Xt £ A,Tq > t), we have to 
write A = {An D) + {AnR*^ X {0}) + {An {0} x W\_ ) and we will compute separately 
the three terms. Let us first remark that Since vi is the unique QSD related to the 
absorbing set dD, it is equal to the Yaglom limit and thus 

F,{X,eAnD) „(Ann) 



^x{TaD > t) 
Let us now write ^1 = A n x {0}. Thus 

F.Mt e ^1) ^ P.(Xt£^i) F,{T^>t,T^=TaD) 
^x{TaD > t) Pi(Ti >t,T^= Tqd) F^{Tqd > t) 

By a similar reasoning as previously, one can prove that p (^t^^^^t'^Tqd) ~^*-*+°° 

z^i,i(^i) and that p (TaD>t) Xil - Ai • ^ similar result inverting indices 

1 and 2 holds for the third term. That concludes the proof of ()4.4p . in the case where 
Ai < Ai^j, for i = 1, 2. 

This result together with the competition case indicates (but we do not have any 
rigorous proof of this) that if Ai = Ai,i for some i, then the QSD (conditioned to 
hitting the origin) is again supported by the axes. 



□ 
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The conclusion of this study is partly intuitive. In the competition case, one of the species 
will kill the other one before dying. In the cooperation case, apparently, all can happen, 
including coexistence of the two types up to a common time of extinction. The problem is 
to know whether all situations for the killing rates are possible or not. 

Here is an heuristic simple argument (to turn it in a rigorous one involves some technicalities) : 
Fix all coefficients equal say to 1 except C12 = C21 = — c with c > 0. Condition (12. 8p reads 
< c < 1. In this situation A12 = An = A and both A and Ai depend continuously on c. 
When c — > (near the independent case) we know that v concentrates on both axes. When 
c — > 1 the rate Ai of return from infinity decreases to (this is the point to be carefully 
checked) so that there is some intermediate value Cc where the phase transition A = Ai 
occurs. 



Appendix A. Killed Kolmogorov diffusion processes and their spectral theory. 

Let D be an open connected subset of M'^, V a C'^ function defined on D. We introduce the 
stochastic differential equation 

dXt = dBt - VV{Xt)dt , Xo = x£D, (A.l) 

for which a pathwise unique solution exists up to an explosion time ^. The law of the process 
starting from x will be denoted by P^;, and for a non-negative measure on D we denote by 

For a subset A of the closure D of D and for e > we introduce 

T% = inf {s > ; d(X^, A) < e} , Ta = lim , (A.2) 

where d{.,.) denotes the usual euclidean distance. T| and Ta are thus stopping times for 
the natural filtration. We shall be mainly interested to the cases when ^ is a subset of the 
boundary dD. 

Our first hypothesis is that the process cannot explode unless it reaches the boundary i.e 

for all X & D, ^ > Tqd , F^. almost surely. (A. 3) 

If D is bounded, ()A.3P is automatically satisfied. If D is not bounded, it is enough to find 
some Lyapunov function. We shall make some assumptions later on implying that a specific 
function is a Lyapunov function, so the discussion on (jA.Sp will be delayed. 

The peculiar aspect of gradient drift diffusion process as (|A.ip is that the generator L defined 
for any function g S C°° by 

L = i A - Vy.V , (A.4) 

is symmetric w.r.t. the measure jj, defined by 

//(dx) = e-2^(^) dx . (A.5) 

The following properties of the process can be proved exactly as Proposition 2.1, Theorem 
2.2 and the discussion at the beginning of Section 3 in [3]. 
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Theorem A.l. Assume that iA.3\} holds. Then there exists a self-adjoint semi-group (Pt) 
on L^(/i) such that for all bounded Borel function f , 

Ptf{x) ■.= ¥.^[f{Xt)lt<Teo]- 
Moreover for any hounded Borel function F defined on Q = C{[0,t], D) it holds 

[F{co) lt<TenH] = Fico) lt<Ten(^) exp (vix) - V{cot) - \ {\VV\^ - /\V){uJs)ds^ 

where E"^^ denotes the expectation w.r.t. the Wiener measure starting from x £ D. 

It follows that for all x £ D and all t > there exists some density r{t,x, .) that verifies 



Ptf{x)=l f{y)r{t,x,y)fi{dy) 
for all bounded Borel function f. 

If in addition there exists some C > such that \ W\'^{y) — AV{y) > —C for all y (z D, then 
for all t > and all x £ D, r{t, x, •) € with 

[ r2(t,2;,y)/x(dy) < (l/27rt)ie^*e2^(^). (A.6) 
Jd 

Remark A. 2. If V{x) +oo as x ^ oo in the condition |VFp(y) — AV{y) > —C for 
all y E D is sufficient for (jA.Sh to hold. Just use Ito's formula with V as in [18] Theorem 
2.2.19. <0> 

We thus have that for any measurable and compact subset A C D and any x D, 

F^iXt £A,TaD>t) = j Fy{Xt-i eA,ToD>t-l) r(l, x, y) fi{dy) (A.7) 

Pt_i(lA)(y)r(l,x,y)/i(dy) 

lA{y) {Pt^ir{l,x,.)){y)fi{dy). 

Hence, the long time behavior of the law of the killed process is completely described by 
Ptril,.,.). 

Since both I^a and r(l, x, .) are in L^(/i), the spectral theory of Pt is particularly relevant. 

This spectral theory can be deduced from the much well known spectral theory for Schrodinger 
operators, thanks to the following standard transform: define for / € Lp{D,dx) (resp. 

Ptif) = e-^'Ptife^) , L/ = 1a/ - l(lVyp-Ay)/. (A.8) 

Pt is then a self-adjoint strongly continuous semi-group on L^(D,dx) whose generator coin- 
cides with L on C^{D). Notice that Pt has a stochastic representation as a Feynman-Kac 
semi-group, i.e. 



Ptf{x) = E^ 



f{Lo{t)) tt<T,o exp (- \ 1^ (iVyp - AV){u;s)ds^ 



(A.c 
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The spectral relationship is clear: if A is some eigenvalue for Pt associated to ip (i.e. Pt{il^) = 
e~^^ij:), it is an eigenvalue of Pt associated with rj = ifj and conversely. 

For simplicity we shall impose conditions ensuring that the spectrum is discrete, in particular 
reduced to the pure point spectrum. Necessary and sufficient conditions for this property to 
hold have been obtained by Maz'ya and Shubin (see [15j ) extending results by Molchanov in 
1953. The criterion is written in terms of Wiener capacity, hence not very easy to directly 
read on the potential V . We shall here assume a less general but more tractable condition 
taken from the Euclidean case explained in chapter 3 of Berezin and Shubin [Ij . To this end 
we now introduce our main hypotheses: 

Definition A. 3. (1) We say that hypothesis (HI) is satisfied if holds and if for all 

Px(raD < +CX)) = 1. 

(2) Hypothesis (H2) holds if 

G{y) = \VV\\y) - lW{y) > -C > -oo 

for all y £ D . 

(3) Hypothesis (H3) holds if 

G{R) = inf {G{y); \y\ > R and y e D} ^ +oo as R ^ oo . 

(4) We say that hypothesis (H4) holds if for all R > one can find an increasing se- 
quence of compact sets Kn{R) such that the boundary of Kn{R) Ci D is smooth and 
lJ^(iCn(i?) n D) = B{0, R) n D, where B{0, R) is the closed Euclidean ball of radius 
R. 

For simplicity we say that (H) holds when and (H2)-(H4) are satisfied. 

We may now state 

Theorem A. 4. Assume that (H) is satisfied, then —L has a purely discrete spectrum < 
Ai < A2 < .... Each associated eigenspace Ei is finite dimensional. If (HI) holds, Ai > 0. 
Furthermore Ei is one dimensional and we may find a (normalized) eigenfunction r]i which 
is everywhere positive. In particular for all f,g in L^(/u), 

hm e^i*(<7,Pt/V = {9,Vi){f,m)- 

t—y+oo 

Proof. The proof of the first statement is similar to the one of Theorem 3.1 in [T], replacing 
B{0,R) therein hy D f] B{0,R). Hypothesis (H4) here is useful to show that the embedding 
B.^{B{0, R)nD) ^ 1?{B{0, R)nD) is compact. Indeed the result is known replacing B{0, R) 
by Kn{R) (due to the smoothness of the boundary). To get the compactness result, it is then 
enough to use a diagonal procedure. 
Ai > is obvious since —L is a non-negative operator. 

It directly follows from the representation formula in Theorem I A . 1 1 ( or (lA.Oh ) that the semi- 
group is positivity improving (i.e. if / > and / 7^ 0, Ptf{x) > for all x G and all 
t > 0). The proof of the second statement (non degeneracy of the ground state ryi) is thus 
similar to the one of Theorem 3.4 in \T\. 

Finally, as in [3] section 3, hypothesis (HI) implies that for / G L^(/u), Ptf goes to in L^(^) 
as t — > +00. This shows that Ai > 0. □ 
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Remark A. 5. Contrary to the one dimensional case, for i > 2 the eigenspaces are not 
necessarily one dimensional. 



Appendix B. Quasi-stationary distributions and Yaglom limit in D. 

The aim of this section is to study the asymptotic behavior of the law of Xt conditioned on 
not reaching the boundary. 

B.l. The general result. The first result is an immediate consequence of the spectral 
theory. 

Proposition B.l. Assume that hypothesis (H) is satisfied. If A C D is such that 1a ^ 
then for all x G D, 

lim e^^* F^{Xt eA,T9D>t) = {1a, Vi) Vi{x) ■ 
In particular if Tji (fi J^^i^ji), limf^+oo P^; {Xt A\T9D > t) = d . 

Proof Recall (|A7l), i.e. ¥^{Xt e A,T9D>t) = f lA{y) {Pt-ir{l,x, .)){y) f^{dy). Since both 
1a and r{l,x, .) are in LP{n) we may apply Theorem I A. 41 and get 

lim e^'^'-^^¥,{XteA,TQD>t) = (1^, r?i) (r(l, x, .), r?i) . 

t— >+oo 

Since r/i is an eigenfunction it holds 

e^i (r(l,x,.),r?i) = e^^Pir]i{x) = rji{x), 

where equalities hold in L^(^). Since r/i satisfies Lr]i = —Xirji in D, standard results in 
p.d.e.'s theory show that r]i is regular (C^) in D, hence these equalities extend to all x & D. 
This yields the first statement. 

For the second statement, choose some increasing sequence D„ of compact subsets of D, such 
that Un Dn = D. It holds 

F.(A. eMTen > t) = l^^^J^'i^'i'^'l < ^^^""^ ^ ^' > 



F^{Xt eD,TaD>t)- P^(At G Dn , Tqd > t) 
so that according to what precedes for all n, 

lim sup ¥.,{Xt £A\T9D>t) < j^^^ ■ 

The infimum over n on the right hand side is equal to as soon as f^^rjidn = +oo, hence 
the result. □ 

In view of what precedes, a non trivial behavior of the conditional law implies that rji G L^(^). 
Conversely this property is enough to get the following theorem whose statement and proof 
are the same as Theorem 5.2 in |3]. Observe that the only thing we have to do is to control 
PtlA for sets A of possible infinite fi mass. 
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Theorem B.2. Assume that hypothesis (H) holds and that rji S L^(/i). 

Then dfi = ijidfj,/ J^r]i{y)fi{dy) is a quasi- stationary distribution, namely for every t > 

and any Borel subset A of D, 

^^,{Xt(^A\TaD>t) = vi{A). 
Also for any x > and any Borel subset A of D, 

^Um^e^^*p,(Ta^ >t) = (^j^viiy)Kdy)^ mi^) , (B.l) 



t—*+co 

This implies that 



lim F^{Xt GA\T9D>t) = ui{A) , 

t—^+CO 



and the probability measure vi is the Yaglom limit distribution. 

Remark B.3. The proof of Theorem 5.2 in |3j Ues on the following estimate 

r{t,x,y) < C{x)e-^''r]i{y) 

for all x,y in D ((0,+oo) in [3]), t > 1 and some function C{x). This result is still true here 
and the proof based on the Harnack's inequality is similar. 



B.2. Ground state estimates. We wish now to give tractable conditions for r]i to be in 
L^(^). Of course if /i is bounded there is nothing to do since rji G L'^(/i), so that this 
subsection is only interesting for unbounded For simplicity of notation we assume that 
the origin € so that ifxED, |x[>a>0. The results of this subsection are adapted 
from Section 4 in [3]. 

Recall that rji = e^ ipi where V'l is the ground state of L (cf. (jA.SP ). i.e. the (positive and 
normalized) eigenfunction of L associated to — Ai. So in order to get some estimates on r/i 
it is enough to get some estimates on ipi . Since V'l = A [tpi ) it is interesting to prove 
contractivity properties for Pi. 

Let us first recall the definition of ultracontractivity. 

Definition B.4. A semi-group of contractions {Qt)t>o "is said to be ultracontractive if Qt 
maps continuously L^(^) in L°°(^) for any t > 0. 

Remark that by duality, and thanks to the symmetry of fj., Qt also maps continuously (fi) 
to L2(^). 

Proposition B.5. Assume that Hypothesis (H2) and liA.9\) are satisfied. Then Pt is ultra- 
contractive. It follows that ^1 = riie~^ is bounded. 

Proof. We may compare the fundamental solution (kernel) of Pt with the one of the Schrodinger 
equation with constant potential as in [Ij by directly using the representation (|A.9p . 
Introduce the Dirichlet heat semi-group in D i.e. 

P,^/(x)=EWM/(a;(t))W,,]= / fiy)p?ix,y)dy. (B.2) 

Jd 
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Hypothesis (H2) and (|A.9p immediately imply that 

where pt denotes the (symmetric) kernel of Pt w.r.t. the Lebesgue measure. 
()B.3P shows that Pt has a bounded kernel for all t > 0, and hence that Pt is ultracontractive. 
In particular ijji is bounded ipi = e^^ sincePi{ipi), hence rji is bounded. More generally 
any eigenfunction ipk is bounded. □ 

From the previous proposition, we deduce that r]i € as soon as /e-^(^)(ix < +00. One 

can improve this result. Recall that denotes the Dirichlet heat kernel defined in ()B.2p . 
Notice that pf{x, y) dy = W^,{TdD > t) goes to zero as x tends to the boundary. 

Proposition B.6. Assume that hypothesis (H) is fulfilled. Assume in addition that there 
exists some R > such that the following is satisfied 

/ e~^^(^')(ix < +00 and / ( / pf{x,y)dy] e"^^"^^ dx < +00. 

J Dn{d{x,dD)>R} J Dn{d{x,dD)<R} \J D J 

(B.4) 

Then r/i € L"^(^). More generally any eigenfunction rj^ € 
Proof. Since rj^ is normalized, 

\ 1/2 

Vk{x)dfi<{l e-^^^^Ux] <+oo. 



/ 



lDr\{d{x,dD)} \J Dr\{d{x,dD)} 

Now 

/ r]kix)dn = / ipk{x)e~^^''^ dx 

J Dn{d{x,dD)<R} J Dn{d{x,dD)<R} 



e^" 7pk{y)pi{x,y)dy] e ^^^^ 



Dr\{d{x,dD)<R} \JD 

< e^/^e^Ml^fc Hoc f ([ p?{x,y)dy) e-'^'^^Ux 

J Dn{d{x,dD)<R} \Jd J 

and the result follows. □ 

Remark B.7. Remark that we can replace by any with s > in the previous proof. 

B.3. Rate of convergence. According to the spectral representation we may decompose 
each function in \? . We introduce some notation. 

Definition B.8. We denote by E2 the eigenspace associated with A2. We know that dim{E2) = 
77,2 < +00 and we may choose an orthonormal basis of E2, {i]2,i, •••,??2,n2)- ^6 denote by pr^ 
the orthogonal projection onto the orthogonal of Mr/i © E2 ■ 

We thus have that 

Pt_ir(l, X, .) = e-^i*r?i(x) + e'^^* ??2,i(x) + e'^^^*"^) h{t, x, .) (B.5) 

1=1 
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where X, .) is orthogonal to Mr/1 ©£'2 and such that II h{t,x,.) ||l2(^)<|| pr-'-r(l, x, .) ||L2(/i)- 
Hence (recall (jXT])) if G ^'^{^J■), 



F,iXteA,TdD>t) = e~^i*(lA,r/i)ryi(x) + e~^2* ^7?2,i(x)(lA, 



V2, 



(B.6) 



1=1 



{h{t,X,.),lA)- 



If we could replace ^ by D we would obtain an expansion of the conditional probability 
^xi^t G A\Tq£) > t). But actually we have 



Lemma B.9. // (H) and are satisfied, Pi is a bounded operator from L°°(/i) to L^(/i). 

Assume firstly the Lemma lB.91 Then we may write 

F^TdD >t)= Pt{lD){x) = Pt_i(Pi(I^))(x) 

with Pi(Iz)) e L2(/i). Note that 

(Pi(Id),%> = / e-^>'r]kdn 
Jd 

since (jB.4p implies that each eigenfunction r/^ is in L^(/i). We thus deduce the 



Proposition B.IO. // (H) and ( B.4 ) are satisfied then for all x ^ D and all measurable 
subset A C D it holds 



t— >+oo 



lim__ e(^^-^i)* (F,{Xt GA\TdD>t)- ui{A)) = 

T.Zim,i{x) ((]iA,r?2,»)(iD,r?i) - {iD,m,i){U,m)) 

r?i(x)((lD,r/i))2 



It remains to prove Lemma lB.91 To this end let us first state an upperbound for pt. 

Lemma B.ll. // hypothesis (H) is fulfilled, there exist a constant M and a non-negative 
function B satisfying lim„^-|-oo B{u) = +00 such that for any x,y in D, 

< pi{x,y) < Me-l^-J'I'/^e-^d^l^lJ'l) . 

Proof. We can obtain an upper bound for pt, when hypothesis (H) is fulfilled. To this end, 
for a non- negative / but this time e = \x\/2 we write 



D 



f{y)pt{x,y)dy 



f{u;{t)) Ii<.^(,) lt<Ten exp (- ^ {\VV\' - AV){us)ds^ 
f{u;{t)) lTen>t>rAe) exp (- ^ {\VV\^ - Ay)K)ds) 



The first term in the sum above is less than 
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For the second term we shall assume that the support of / is included in the ball B{x,e/2). 
Recall that for a brownian motion starting at x, the exit distribution from B(x, e) is uniform 
on the sphere S{x,e). Hence 



S{x,e) J 



dy 



2 \ -d/l 



where 7 is the ordinary heat kernel. Since I2; — y| > e/2 in the above formula, 

< (2vru)-'^/2g-^V8«<g-d/2 f^^^^ 

the latter inequality being obtained by an easy optimization in u. Since |x| > a this quantity 
is bounded on D by some constant B. But 

for some constant K depending on d only. So 

and gathering all the previous results we obtain 

if < |x|/4;pi(x,y) < M ((2vrt)-'^/2 e'^^d^l/^V^ + g-l^l'/^^*'^) , (B.7) 

for some constant M . If — y| > |x|/4, (IB.3j) furnishes 

Pt{x^y)< e^*/2(27rt)-'^/2e-l^l'/32*. 

Define 

B{u) = ^ min (G(n/2)/2 , {u^ /?,2d)) . (B.8) 

We have shown that there exists some constant M such that pi{x,y) < Me^^^^l^'^ but 
since pi is symmetric the same holds replacing \x\ by \y\ and finally |x| by max(|x|, \y\). 
Taking the geometric average of this estimate and ()B.3p ends the proof. □ 

Proof, of Lemma IB. 91 let us consider a bounded function H^Hoo < 1- Then Pig{x) = 
e^(^) Pi(e~^5f)(x). According to ()B.4p and lemma [BlT| 



/ \Pi{e-^ g){x)\dx < 

J Dn{d{x,dD)<R} 
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J Dr\{d{x,dD)<R} \JDr\{d{y,dD)>R} J 



+ M I If e-^^y^p^{x,y)dy] dx 

J Dn{d{x,dD)<R} \J Dn{d{y,dD)<R} J 



/ \ 1/2 

< M I [ e^2y(,)^y 

\JDn{d{y,dD)>R} J 



Dr\{d{y,dD)>R} \J Dr\{d{x,dD)<R} 



.2 \ 1/2 
g-|x-yp/4g-B(|x|vM)^^] 



+ M [ e-^(^) ( [ pf{x,y)dx) dy 

J Dn\d(y,dD)<R\ \Jd J 



' Dr\{d{y,dD)<R} 

is finite (recall that y) =Pi{y,x)). Hence l^^^^ Q£)-^<jiPi{e~^ g) ^L^{dx). Since A is ul- 

tracontractive, Pi is a bounded map from ^J-{dx) to lJ'{dx). It follows that Pi{'^d[x,dD)<R Pi{^~^ 9)) ^ 
L'^{dx). 
In addition 

! \Pi{e~^9){x)\^dx< f e-2^(^)dx < +00 , 

JDn{d{x,dD)>R} J Dr\{d(x,dD)>R} 

i-e. Id(x,5r>)>i? A(e-^fif) e ^^(c^a;) that A(Id(x,OD)>i? A(e"^5')) e ^^(da;). 
Summing up yields that P2{e~^ g) G ^^{dx). 

Actually we may replace 1 by s > as remarked in remark \Q.7\ thus replace 2 by 1 in the 
previous result. □ 



B.4. Uniqueness of the quasi-stationary distribution. In [3] Theorem 7.2 we derived 
a necessary and sufficient condition for ui to be the only quasi-limiting distribution, i.e. to 
satisfy limf_^+oo IPf (^t ^ ^\Tqd > t) = I'liA) for all initial distribution i'. In that case vi is 
the only quasi-stationary distribution. (Recall that in D = This condition is very 

close to the ultracontractivity of the semi-group Pt. 

We shall not try here to obtain such a criterion, but only a sufficient condition based on the 
previous remark. 



Proposition B.12. Assume that (H) and {B.4) are satisfied. If Pt is an ultracontractive 



semi-group, for all initial distribution v and all Borel subset A C D 

lim F^{Xt eA\TdD>t) = iyi{A) . 

t—*+co 

In particular vi is the unique quasi- stationary distribution. 



Proof. Since Pt is ultracontractive, it turns out, as proved in [6] Theorem 1.4.1, that ^^(/u) fl 
L°°{ijl) which is included into L^(/i), is invariant under Pt- So Pt extends as a contraction 
semi-group on all L^(^). 
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Now, according to Theorem I A . 1 1 we know that ¥iy(Xt G A , Tqd > t) = r^{t,y)ii(dy) with 
r^{t, ■) = f r{t, X, .)h'{dx) £ and r^{t + s,y) = Ps{ru{t, ■)){y). Hence for t > 2, 

F,{Xt e A , Ton > t) = [ lAPt-iir,{l,.))dfi 

Jd 

= [ Pi{lA)Pt-2{r,{l,.))dfi. 
Jd 

But thanks to Lemma lB.91 and to ultracontractivity, both Pi(I^) and Pt^2iru{^, ■)) are in 
L^(/i). Furthermore e"^^^^^^^ Pt-2ifu{^, ■)) converges strongly in to (Pi(r,y(l, .)), 771) r/i as 
t — > +00. Hence 



hm ¥,{Xt eA\TeD>t)= p / ' " = MA) 



since rji G L-'^(^). 



□ 



It remains to give tractable conditions for Pt to be ultracontractive. To this end we shall use 
the ideas introduced in [13j and later developed in [2] in particular. The next lemma is the 
key 



Lemma B.13. (see |13j ) Assume that (H) and {B.4) are satisfied. If for all t > there 
exists c{t) such that for all x £ D, 



< cit) 



(B.9) 



then Pt is ultracontractive. (Recall that G{y) = |Vyp(y) — lSV{y)). 

Conversely this condition is necessary if we assume in addition that e^^^^ fi{dx) < +00. 

Proof. The proof is the same as in It is given for the sake of completeness. 

Recall that the heat semi-group on D is ultracontractive, i.e. for all non- negative / G L^((ix) 

supE^- [lT,^ytf{uJt)] < (vrt)-^ II / ||L2(rf.) . 
x&D 

For a non- negative g G L^((i/i), / = e^ g G lJ^{dx) so that using Theorem lA.il 
P,g(x) = e^(-)EW^ ^' . -i r*, 

< e^(^t)-i \\g\k2^,,) e^(-)E^ 
Hence if ()B.9P is satisfied, 

Ct 1 

Pt9{x) < c{t/2)e- (7rt)-4 || g 

for all X £ D, i.e. Ptg is bounded and Pt is ultracontractive. (jB.OP is thus a sufficient 
condition. It is also necessary once e^ G L^(/i), since 



Irs^>i/(^i)e-HoGKK 

-TBD>t/2 



t/2/K/2ie 2 JO 



Pt{e^){x) = e^(^)E^ 



lTsD>te 



-5 Io'^i^s)ds 



as it can be observed in the Girsanov formula stated in Theorem lA.il 



□ 



26 P. CATTIAUX AND S. MELEARD 

Papers [13] and [2] contain several methods to prove (|B.9p . The most adapted one to our 
situation is the " well method" based on the Girsanov transform (Theorem lA.ip . 

To this end we shall introduce some notation : 

• ior R>0, AR = Dn {\x\ < R}, 

• V{R) = supj-g^^ V{x) (be careful that there is no absolute value), 

• for e > 0, De = {y e D , d{y, dD) > e] and Se = inf{t > , uj{t) e Df], 

• for /c G N* , ek = inf{i > , uj{t) e Ak}. 

We then have the following analogue of Theorem 3.3 in [13j 

Proposition B.14. Assume that (H2) and (H3) are satisfied. We shall also assume 

that for all e > {),V is hounded from below on Z?^. Let a^. = G{k), = V{k) and 7^ be such 
that YlT=i 7fc < +00. Pt is ultracontractive as soon as the following holds: 

/ora///?>0, ^expQ (6fe+i-/37fcafc)^ < +00 . (B.IO) 

Proof. The first step is to be convinced that Lemma 3.1 in [13j is still true i.e if r/j = inf{t > 
, G Ar} and x ^ Ar , 

^vi.) JEW. hr,^^,^ e-5 /o^^ < ^v{R) 

Define Mt = e^^^"^*)-^ IoG{^s)ds_ Thanks to our hypothesis on V and (H2), MtAmASe is 
actually a bounded martingale. The result follows by making successively t go to infinity and 
e go to 0. 

Once this is proved the rest of the proof is exactly the same as in [I3j except that we have 
to replace the stopping times therein by A Ss and then make e go to again. □ 
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